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1 Introduction 



The process of diffractive p meson electroproduction 

7>->PP (1.1) 

is now under intensive experimental study at the ep collider HERA in a new 
energy range as compared to the previous fixed target experiments. The word 
diffractive means that the energy of ^*p collision W is much larger than the 
virtuality of the photon Q. Therefore this process is mediated by the pomeron 
exchange. On the other hand the photon virtuality is large in comparison with 
the typical hadron scale and we deal here with hard diffraction. 
The diagram of this process and our notations are given at Fig.l 




Figure 1: Diffractive electroproduction of vector meson on a proton target. 

Theoretically this process was considered both in nonperturbative and per- 
turbative approaches in many papers. The references can be found in the review 
PHI . It turns out that the perturbative models where pomeron is represented 



by the hard two-gluon exchange are able to reproduce the main features of the 
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HERA data for hard diffraction. In these models the amplitude is proportional 
to the gluon distribution in the proton 

M oc xG(x,Q 2 ) , (1.2) 

where x = 0%, Q is the typical hard scale of the process. 

The aim of this paper is to consider in the framework of perturbative QCD the 
polarization effects in diffractive electroproduction of vector meson at large Q 2 . 
Experimentally, information about the polarization state of produced meson is 
extracted from the angular distributions of the meson decay products (7r + 7r~ for 
p). The present day results of such investigations at HERA are consistent with 
the s- channel helicity conservation (SCHC), which means that the produced 
meson retains the helicity of incoming virtual photon. The expected increase 
of the HERA luminocity gives us a hope that we will have in the near future 
much more precise experimental information about helicity properties of the 7*M 
transition. Such an information will give us valuable insight into the underlying 
dynamics. 

We shall discuss experiments with unpolarized protons, therefore the proton 
can be formally considered as a spinless target and we shall indicate further only 
the polarization states of the virtual photon and the produced meson. 

Under the asumption of SCHC there are only two independent helicity am- 
plitudes M 1 * l ^ Pl and M 7 *^ PT , where L(T) denotes longitudinal (transverse) po- 
larized state. 

The longitudinal amplitude dominates at large Q 2 . It has been considered 
in the papers fl], ||, ||, [§, ||. A QCD factorization theorem for this amplitude 
was proven in ||. The situation for the transverse amplitude is more compli- 
cated. Formal power counting gives ~ ^ suppression in comparison with the 
longitudinal amplitude. However, this suppression factor is too small for typical 
values of Q 2 ~ 10 GeV 2 for HERA electroproduction experiments and does not 
agree with the measured ratio of For the transverse amplitude the integral 
over the longitudinal fraction z of the quark momentum is logarithmically di- 
vergent in the end points of the integration region. This means that essential 
transverse distancies which scale like — , 1 can be large even for large Q 2 and 

that the transverse amplitude can receive large or even dominant contribution 
from the nonperturbative region. As was noticed in 0, if the nonperturbative 
contribution would be dominant we would expect for cry features similar to the 
ones observed in photoproduction experiments as: a "softer" as compared to 
energy behaviour ~ W ' 2 and a larger slope b ~ 9 — 10 GeV -2 . However these 
expectations are not supported by the data. 

In || it was assumed that the scale behaviour of gluon distribution G(x, Q 2 ) ~ 
(qz) , where 7 is the anomalous dimension of gluon density, plays a very impor- 
tant role in the physics underlying the transverse amplitude. Taking into account 
this dependence it can be seen that the typical transverse distances are ~ 
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and can be smaller than A( ^ CD at large Q even for small 7. Therefore perturbative 
QCD can be applicable to the transverse amplitude. The estimate derived in || 
under the assumption of constant anomalous dimension 



a T M p 2 l v l + 7/ 



1.3) 



shows the role of the scaling violation of the gluon density and agrees qualitatively 
with the data. 

Under the assumption of natural parity exchange in the t— channel there 
are five independent helicity amplitudes. These are the two helicity conserving 
amplitudes 

M (0 ,o) = M % _ po , M (+li+1) = M^ p+1 (M(-i-i) = M (+li+1) ) ; (1.4) 

and the amplitudes violating SCHC: two single spin-flip amplitudes 

M (+1)0) = M lU ^ po , (M ( _ lj0) = -M (+1)0) ) , (1.5) 

M (0 , +1) = M^ p+1 , (M( _!) = -M ( o,+i)) ; (1.6) 
and one helicity double-flip amplitude 

Mm-v = M r+i ^ , (M ( _ li+1) = M (+w) ) . (1.7) 

As it will be shown below, similar to the non-flip amplitudes ( |1.4|) , the single 
spin-flip amplitudes can be expressed through the gluon density in the proton. 
In the leading order of 1/Q expansion the double spin-flip amplitude does not 
receive a logariphmic contribution from the integration over the t— channel gluon 
momenta, see Fig. 1. Therefore in the first term of the 1/Q expansion of the 
double spin-flip amplitude the large factor xG(x, Q 2 ) is absent. 

We will show that the largest amplitude violating SCHC is M( + i i o). The other 
amplitudes violating SCHC can be neglected in a first approximation. Our result 
derived in the approximation of constant gluon anomalous dimension 7 reads 



(3=^7^ = ^— ■ (1.8) 



M( ,o) V2Q 
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The observed t— dependence for p meson electroproduction is ~ e , with 
the slope b = 5 ... 6 GeV~ 2 . At typical values of t ~ 1/6 this amplitude is not 
too small and, as will be shown below, leads to a sizable interference effect. 

In the perturbative QCD approach the effect of the scale behaviour of the 
gluon density manifests itself qualitatively in the similar way in both transverse 
amplitudes: the helicity non-flip M(+i j+ i) and the helicity single-flip amplitude 
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M( + i 5 o) . Therefore the mesurement of Mr+i t o) at HERA would give us an impor- 
tant check of whether perturbative QCD describes correctly the physics under- 
lying the amplitudes of vector meson electroproduction initiated by a transverse 
photon. 

The present work is based on the experience gained in a previous study of the 
diffractive vector meson production 0. 

The paper is organized as follows. In section 2 we present the main steps of 
our calculation method. The results for helicity amplitudes and the discussion 
of the underlying physical effects are given in section 3. The influence of the 
helicity-flip amplitudes on the vector meson product distributions is discussed in 
the section 4. Our conclusions are summarized in section 5. 

2 Impact parameter representation and the me- 
son wave functions 

The amplitude of the diffractive process 7*p — > pp can be represented as the 
integral over the transverse momenta of gluons in the t— channel (impact repre- 
sentation) 

d 2 k 



M 7 * p ^.pp — iSry* p J ^ 2 ^q _ ^^J'Y'pJp ■ (2-1) 

Here q is the momentum transfer which is transverse with high accuracy, q 2 = t = 
— q 2 . Throughout the paper all vectors, if it is not mentioned separately, are two- 
dimensional vectors in the transverse space. The accuracy of the representation 
(|2.1|) is expected to be ~ Q 2 /s 1 * p . 

The space-time picture of the process in diffractive (high energy) region is 
the following. The virtual photon fluctuates into the qq pair long time before and 
the qq pair converts into the vector meson long time after the interaction with 
the proton. Therefore it is possible to represent the photon impact factor as the 
convolution of the impact factor for the qq dipole scattering with the light cone 
wave functions of the incoming virtual photon and the outgoing vector meson 
(Fig.2) 

■Wk,q) = / <Mli, *i)^°' e (li, la, *i, to, k, q)%(h, z 2 ) , 

(2.2) 

where 

Qdipole^ ^ ^ ^ q ) = 167r 3sfl^±^ JS _ Zi) _ lx _ qz ) + 

5(1 -h + qz) - 5(\ - h + k - qz) - 5(\ - h - k + qz)} . (2.3) 

iV = 3 is the number of colors, (5 ab ) 2 = N 2 — 1. as is the strong coupling 
constant. 
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q-k 



Figure 2: The 7* — > M impact factor. Only one contribution of the interaction 
of the qq dipole with the two exchanged gluons (dashed lines) is shown. 



The light cone wave function of the photon 



uev 



-e q \j zz- 



'\ 2 + Q 2 zz (2 ' 4) 

describes the probability amplitude for the splitting of the photon into the qq 
pair with electric charge e q , (e q = |e for the u quark). The quark carries the 
transverse momentum 1 relative to the photon momentum and the fraction z of 
the photon longitudinal momentum (the fraction for antiquark is z — 1 — z). 
It should be noted that the vector 1 2 in eq. ( j2.2|) is transverse relative to the 
momentum of the outgoing vector meson which received non-zero momentum 
transver q. 

The polarization state of the photon is described by the vector e. We choose 
the following convention for the polarization four vectors describing a transversly 
polarized vector meson and a transverse photon 

e±=^(0,l,±i,0), (2.5) 

where the x axis is choosen in the direction of the momentum transfer 

q = g- (0,1,0,0). (2.6) 

It can be seen that the difference in the polarization vectors of the vector meson 
from (|2.5| ) related to the non-zero value of momentum transfer can be neglected. 
The polarization vector of the longitudinally polarized virtual photon is 

e = ^(pi,0,0, v /p?-g 2 ) , (2.7) 

where p± is the value of the longitudinal momentum of the photon. And the 
similar convention is adopted for the polarization vector of the longitudinaly 
polarized vector meson. 
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The perturbative expression ( |2.4| ) for the splitting of a photon into the light qq 
pair has the important property that the total helicity of the produced massless 
quark-ant iquark pair is zero irrespective of the polarization state of the photon: 
a) for longitudinally polarization 



'zzu\eoV\> = 2Qzz5 x ,~x' , (2.8) 
b) for transverse polaritation 

sfzlu x e±v x , = 5 A) _A'{(1 " 2z) =F A}(e±l) , (2.9) 

The quark helicities ±| are represented in the above equations by A = ±1. 

The splitting of the transverse photon into the quark pair with the total 
helicity ±1 is proportional to the current quark mass. Since light quarks have 
very small current masses we will neglect here this splitting. Of course it becomes 
important in the case of a heavy flavour production. 

The physics of the spin flip transitions which is the main subject of this paper 
looks more transparent in the space of impact parameters of the qq pair. The 
helicities of quarks coincide in this case with the projection of the quark spins 
onto the z axes. Since the total helicity of the quark pair is zero, the projection 
of the orbital momentum of the incoming pair onto the z axis should coincide 
with the helicity of the incoming virtual photon. Another important property of 
perturbative QCD is that the interaction of t— channel gluons with the pair does 
not change with high accuracy the helicity states of the quarks. This interaction 
does not change also the impact parameters of the pair even if the momentum 
transfer is not zero. Therefore the helicity state of the produced meson should 
coincide with the projection of the orbital momentum of the outgoing quark pair 
onto the z axis. We can conclude that in the frame of perturbative QCD the 
only posibility to have the change of the helicity state during the diffractive 7*p 
transition is to change the z projection of the angular momentum of the qq pair 
in the interaction with the proton. 

To make this discussion more quantitative let us transform eqs. (|2.2| . |2~"3 . 



2~^) into the space of impact parameters. By Fourier transformation we obtain 



the photon (meson) wave function in the representation of impact parameters 

*(t,z)= [ p-9(l,z) - e- ilr , (2.10) 
where r is the difference between the quark and antiquark impact parameters. 



^.(r, z) = TieQy/zz 



Vv{(l - 2z) =F \}{e±-)K x {rQ^zz) 

r 



(2.11) 



Vt(r,z) = -2eQzz5 X) - X ,K a (rQ\/zl 



(2.12) 
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In these eqs. ( |2.11|J2~T2| ) Ko : i(mr) are the Bessel (MacDonald) functions. 



We can represent the impact factor describing the 7* — > p transition (2.10) in 
the form (Fig. 2) 

■Wk,q) = / ^^* r (r 1 ^)^(r 1 ,r a) k,q,z 1 ,^)*;(r 2> z)^^ 

^(2.13) 

where $ d *P° Ze (r l5 r 2 , k, q, Zi, z 2 ) is the Fourier transform of ( f2.3|) 



$«( ri , r, k, q, Zl , z 2 ) = 167r 3 ^5(r - n)<^i - z 2 )f{k, q, r, Zl ) , 

/(r, k, q z) = e iqrz (l - e~ ikr )(l - e^( q - k ) r ) . (2.14) 

The expression for the dipole impact factor is proportional to S(r — r x ) which 
reflects the mentioned above property that the interaction does not change the 
impact parameters of the pair. The main part of the dipole impact factor is the 
factor /(r, k, q, z). This factor tends to zero if the momentum of one of the t— 
channel gluon (k or q — k) or if the transverse separation between the quarks r 
vanishes. 

Let us discuss now the wave functions of a vector meson. They can be con- 
structed using the analogy with the photon wave functions. The wave function 
of the longitudinally polarized photon can be rewritten in the following way 

Where M is the invariant mass of the qq pair. We shall adopt the natural as- 
sumption that the meson wave functions depend on the invariant mass of qq pair. 
But in the p meson case the corresponding $(M) should fall off faster at large 
M as compared to the photon wave function. Therefore 

*,M=- 3 ^t^&®-. (2.16) 

The dimensionful coupling constant f p ~ 200 GeV is related with the e + e~ decay 
width of the p meson, 

T = J -£- , a = e 2 /4n = 1/137 . (2.17) 

3m p 

The normalization for function $(M) is 

d 2 l I 2 

:$( — ) = 1 . (2.18) 



(16tt 3 )zz zz' 
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The wave function of the transversly polarized p meson 



zz{(l - 2z) =f A}(e± 



l.$(l 2 /(zz)) 



zz 



(2.19) 



is constructed in analogy with the wave function of the transverse photon. In 
the meson rest frame the meson splitting looks like as follows: the quark and the 
antiquark fly out back to back, the polar angles of the quark momentum relative 
to the direction of the meson momentum in the boosted frame are (#,</>). The 
amplitude of the transition of a spin 1 p-meson in a state with the helicity A 
into the qq pair with the total helicity A' is proportional to the rotation matrix 
d\ x ,(8, 0). We are interested in such meson splittings when the quark and anti- 
quark have the opposite helicities (zero total helicity) in the boosted frame, that 
corresponds to the total helicity of the pair ±1 in the rest frame. Since the angle 
9 in the rest frame and the variable z in the boosted frame (the fraction of meson 
momentum carrying by the quark) are related as 

1 + cos(fl) 



it can be seen that the ratio of ^ (1, z) and ^(1, z) is equal to the ratio of the 
corresponding rotation matrices. This observation justifies eq. ( |2.19| ) for the 
wave functions of transverse p as well as the relative sign between the transverse 
and longitudinal wave functions. 

The meson wave functions in the representation of the impact parameters are 
given by the Fourier transform of eqs. ( |2.16| , |2TT9| ). Since the qq fluctuation of the 
transverse polarized meson has the projection ±1 of the angular moment onto 
the z axis its wave function has a factor (e±r). 

Now let us return to eq. (|2.13|) and discuss the virtual photon to meson 
impact factor 



J r _ p (k,q) 



a s S ab rd 2 rdz 



N J 16tt 3 j 
for the various helicity transitions. 



# 7 *(r, z)/(r,k, q, z)**(r, z) 



(2.20) 



Looking at the eqs. (|2 . 1 6| , pTT9| , P . 1 1 | , pTT2|) for the meson and photon wave 
functions it is easy to see that: 

a) the transitions obeying SCHC are proportional to the dipole factor /(r, k, q, z) 
averaged over the polar angle of the vector r 



(/ (r, k, q, z)) = [J {rqz) - J (r\k -qz\)] + [z 



<-> z\ 



(2.21) 



b) the single spin-flip transitions are proportional to the projection of this dipole 

k — qz 



factor onto £ 



(-/(r,k,q,z)) 

r 



-Jx(rqz) 
Q 



qz\ 



- Ji(r|k — qz\ 



[z «-> z) 



(2.22) 
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c) the double spin-flip transitions are related to the projections of this dipole 



factor onto 



r u iv /r u r 



-f(r, k, q, z)). They are proportional to 



2 (k q^q) 2 j2(r|k _ j} _ q 2 (J (r j k _ j } + Urqz)) 

(k — qz) z 



[z^z] . (2.23) 



Using the general formulae derived in this section we shall calculate in the 
next section the helicity amplitudes in the high Q 2 limit using the approximations 
to be described in the following. 



3 Calculating the helicity amplitudes 

Assuming some functional form for the meson wave function allows to calculate 
according to the eq. (|2.20 ) the 7* — > p impact factor without any approximation. 



Unfortunately the meson wave functions are purely known. But, as it will be 
seen further, if we proceed to calculate the amplitude in the leading order of 4 
expansion only a limited information about these wave functions is needed. 
The eqs. (|2.11J2.12| ) show that the typical size of the virtual photon fluctua- 



tion is r ~ Q^Tji • Outside of the end point regions of z this size is much smaller 
than the meson transverse size. Therefore to calculate the leading in ^ behaviour 
of the impact factor we have to calculate the first term of the Taylor expansion 
for the meson wave function and the expansion of the dipole factor /(r, k, q, z) 
in the region of small r. 

Making the Fourier transform of eqs. (|2.16|JQ9l) and then expanding them 
at small r we find 



*J(r, z) w -^—"-5^f p zz , (3.1) 



16tH3 



p 



167T 3 3 



2tt 2 



where 



* (r,2)«TiT" -5 x ,-wf P zz{{l-2z)^X\{e ± r){M) , (3.2) 

is the mean invariant mass of the qq fluctuation. This mean invariant mass is 
expected to be of order of the p meson mass. We would like to note that the 
Taylor expansion for the transversly polarized meson wave function starts from 
the term ~ r, whereas the longitudinal polarization wave function is constant at 
small r. The appearance of this suppression factor r is related with the nonzero 
projecton of the orbital momentum (±1) of the quark pair in the transverse case. 

The azimuthal projections of the dipole factor /(r, k, q, z) for the various 
helicity transitions , see eqs. ( |2.21| - [2~2"3| ), can be simplified in the region of small 



r by expanding the Bessel functions. 
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a) for the SCHC transitions: 

4 



Wx) » 1 - €) 



(f(v, k, q, z)) = ^(k 2 - (kq)) + 0(r 4 ) ; (3.4) 



b) for the single spin-flip transitions: 
(Ji(x) « f - 

£/(r, k, q, *)) = » r3(g 2 7 z) [k(g 2 - 2(kq)) - q(A; 2 - 2(kq))] + 0(r 5 ) ; (3.5) 

c) for the double spin-flip transitions: 

-dj [2(kq) 2 - (kq)g 2 - k 2 q 2 } 

[2A; 2 (kq) 2 - jfcV - 2(kq) 3 + 



4 • 12g z 

3(kq) 2 gV + z 2 ) - 2(kq)g 4 (z 3 + 2 3 )] + 0(r 6 ) . (3.6) 

We quote in the last case the expansion up to the next-to-leading term, the 
importance of which will be discussed in what follows. 

Let us discuss now the proton impact factor and the integration over the t- 
channel gluon momenta in the eq. ( |2 . 1| ) . 

The proton is a colorless state. Therefore its impact factor vanishes if the 
transverse momentum of any of the t- channel gluon tends to zero. On the other 
hand if the transverse momenta of gluons are large, much larger than the inverse 
transverse size of the proton (the value of the momentum transfer q is expected 
to be small, therefore k w k — q in this region), both of the gluons couple to 
the same parton inside the proton and as a function of k the impact factor is 



approximately a constant in this region. According to eqs. (3~3,^5") the impact 
factors of the helicity conserving and single helicity flip 7* — > p transitions are 
proportional to the square of the t- channel gluon momentum k 2 at large k: 

a) for the SCHC transitions 

</(r,k,q,z)>^^; (3.7) 

b) for the single spin-flip transitions 

C-firX^z))*- ^ 3 ^*-^ . (3.8) 

Therefore the main, logarithmic, contribution to the helicity non-flip and helicity 
single-flip amplitudes originates from the broad region of large k, k < Q = K^ . 
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In the leading log Q 2 approximation and at t = these amplitude are proportional 
to the gluon distribution 



o 2 

ftab 7 J 1.2 
\2\ / T ab a ^ 



x-G(x,Q') = -J J?^ . (3.9) 
Let us calculate these amplitudes in the leading log approximation. 
3.1 SCHC amplitudes 

Let us start with the dominant at large Q 2 longitudinal amplitude. Using ( |3.9| ) 
and inserting eqs. (|2.12| , |3~T|J3.71 ) into eq. Q2.20| ) and performing the sum over the 



quark helicities we obtain the following expession 



M( 0)0 ) = is 7 * p 



J d 2 rdzr 2 3 ^ S J p Q(zz) 2 K (rQ^)xG(x,Q 2 zz) . (3.10) 



We took into account that the mean electric charge of the quarks inside the p 
meson (|p) = -j^(\ u u) — \dd))) is e/y2. Performing the integral over r, using 

oo 

/ Ko(r)r 3 dr = 4, we obtain the result 
o 

M (0j o) = i'Vp / dz & ^nq3^ xG ( x, Q 2z ^ ■ ( 3,n ) 

If we neglect the z dependence of the argument of the gluon density and take the 
integral over z in eq. (3.11) we will reproduce the known result for the longitudinal 



amplitude QTJ, ( for the asymptotical form of the meson distribution amplitude). 
The integral over z is convergent in eq. ( |3.11| ) . Therefore the end point regions 



of z do not bring the essential contributions, and typical transverse distances of 
the process are small, ~ ( q^= )- This justifies the application of the perturbative 
QCD in this case. 

The situation is more complicated for the case of the transverse amplitude. 
Inserting eqs. ( CTJE^JTTD into Q ) we obtain 

M i+h+1) = iSrp Jd 2 rdzr 3 ^^ 
(M)Q{zzf/ 2 {z 2 + z 2 )K 1 (rQ^)xG(x, Q 2 zz) . (3.12) 



The integration over r gives, using / Ki(r)r A dr = 16, 

o 

M(+i,+i) = is rp J dz ^2 N q 4 ^ ( z + z )xG{x,Q zz) . (3.13) 
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If we would neglect the scale dependence of the gluon density G(x, Q 2 = Q 2 zz) the 
integration over z would be logarithmically divergent at z — > 0, 1 in the above 
integral. On the other hand in the small x region the gluon density increases 
rapidly with Q 2 . This increase can be simulated in a first approximation as 



G(x,Q 2 zz) 



G(x,Ql) 



(Q 2 zz)/Q 2 



(3.14) 



with a constant anomalous dimension 7 of the gluon density. Then the integral 
over z in the eq. ( |3.13| ) is convergent and we confirm the result of || 



a 



M (+ i, + i) _ (M) 1 + 7 



Mi 



(0,0) 



Q 



7 



(3.15) 



3.2 The helicity single— flip amplitudes 

The calculations of the two independent helicity single-flip amplitudes are quite 
similar to the ones for the SHCH amplitudes. 

Using eqs. ( |2.11| , |3~T| , |3.8| ) we obtain the following expression for the single 
spin-flip transition of the transversly polarized initial photon 



M [+i,o) = is rp J d 2 rdzr 3 ^^ 
\t\Q(zzfl 2 (z - z) 2 K 1 (rQy/zl)xG{x, Q 2 zz) . (3.16) 
The integration over r gives 



r 2 ■ 37c 2 easfpJ\t\ 
M(+i, ) = i'Vp J dz NQ^^zz) — ^ ~ Z > xG ( x iQ zz ^> ■ ( 3 - 17 ) 

According to eqs. ( |2.12| , |3.2| , |3.8|) the single spin-flip transition of the longitu- 
dinaly polarized initial photon is 

M (0 ,+i) = -is rp J d 2 rdzr 4 



\t\{M)Q(zz) 2 (z - z) 2 K (rQy/z^)xG(x,Q 2 zz) . (3.18) 

00 

After the integration over r, using J K (r)r 5 dr = 64, we have 







f 4-3n 2 ea s fJ\t\{M) 
Af(o,+i) = -«Vp J dz NQ 5 (zz) ^ ~ ^ xG ( x, Q ZZ > ■ ( 3 - 19 ^ 

We calculate the integrals over z in (|3.11| , |3.17| , |3.19|) using the assumption of 
the constant gluon anomalous dimension ( |3.14| ). The following results for the 
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ratios of the single-flip amplitudes to the longitudinal non-flip amplitude can be 
derived 



P 



M, 



(+1,0) 



Mi 



(0,0) 



Q V2-f 



Mi 



(o,+i) 



(M)J\t\y/2 



Mi 



(0,0) 



Q 2 
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(3.20) 
(3.21) 



3.3 The helicity double— flip amplitude 



The helicity double spin-flip amplitude Mt + \ -\\ can be written as the sum of the 
two parts 

^(+i,-i) = M(° + i,-i) + M\ +Xt _ x) , (3.22) 
where M? +1 _ 1 n corresponds to the first 

corresponds to the second (~ r 4 ) term of eq. ( |3.6|) . Keep in mind that the 
expansion in r is equivalent to the expansion of the amplitude in 



'~ r 2 ) term of eq. ( |3.6| ), and M} +1 _ x s 



The integrations over the transverse momenta of the t— channel gluons, see 
eq. (2TI), and over the qq longitudinal momentum fraction z are diffrerent for 
Mf+L-U and Ml +lt _ iy 

Let us discuss M? +1 _^ first. Using ( |2.11j , |3.2|,|3.6|) we obtain 



M? +1 



-is rP J d2r ^f ^ f ^ {M)Q{zzfl 2 K 1 {rQVrz)h , (3.23) 



where the factor I represents the integration over the transverse momenta of the 
t— channel gluons 



a s 5 



ah 



d 2 k 



q 2 J k 2 (k — q) 



2(kq)^ - *Y - (kq)g^ J" fe (k,q) . 



(3.24) 



Performing the integration over r in (|3.23| ) we have 



M (°+i,-i) 



i s 



TP 



V2NQ 4 



(3.25) 



The above equation should be compared with the corresponding eq. (|3.13|) for 
the helicity non-flip transverse amplitude M( +li+1 ). In contrast to ( |3.13 ), the 
integration over z in ( p.25|) is not singular in the end point regions. This difference 
is related with the helicity flip. Looking at the expressions ( |2.11| ) and ( |3.2| ) 
for the photon and vector meson wave functions it is seen that in the flip case 
the sum over the helicities of the intermediate qq pair is proportional to the 
additional factor zz. For the non-flip case the sum over qq helicities gives the 
factor ~ (z 2 + z 2 ), which does not vanish at z = 0, 1. 
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Now let us discuss the integration over the t— channel gluon momenta ( |3.24j ). 
This integral is convergent on the upper limit. Therefore M? +1 _^ can not be 
expressed through the gluon density using eq. (|3.9|) . Moreover, since the main 
part of the integral ( |3.24|) originates from the region of small transverse momenta 
of the t— channel gluons, k ~ q, we are dealing here with soft physics. 

We shall assume that this soft t— channel exchange can be described in the 
frame of the two-gluon exchange model, see eq. (|2.1| ), with some simple functional 



form for the nonperturbative impact factor of the proton J p . The model of such 
type was used successfully in [|J to describe the total pp cross sections. Following 
0, the proton impact factor is 



J p (k,q) = a s S 



ab 



A 2 



A 2 



A 2 + q 2 /A A 2 + (k - q/2) 2 



(3.26) 



The first term in ( |3.26| ) describes the contribution of the diagrams where both t— 
channel gluons are coupled to the same quark inside the proton. This contribution 
is similar to the electromagnetic form factor, therefore it is natural to adopt that 



.4 



m. 



(3.27) 



The strength of the nonperturbative coupling a>s = f- is a free parameter. To 
choose its value, let us calculate the value of the total pp (pp) cross section which 
is related through the optical theorem to the forward amplitude 



pp 

a tot 



Im{MPP(t = 0)) 



(3.28) 



3 pp 



Using eqs. (|2.1| ) and (|3.24 ) it can be shown that 



Im(M pp (t = 0)) 



A 2 



(3.29) 



Describing the pp total cross section with this simple model leads to an parameter 
as increasing with energy. We understand that one can do better by replacing 
the two gluon exchange by a pomeron exchange. For our aim of a simple estimate 
we shall adopt the relations ( |3.28| , |3.2^ ). 

Performing the integral ( |3.24| ) (using the expression ( |3.26| ) for the proton 
impact factor) we find 



A 2 



A 2 + q 2 /A 



A A 2 n 2 

:i log (i +-£--) 



4A 2 ' 



(3.30) 



At small q 2 , q 2 « A 2 



71 



4|t| 



A 2 



7T- 



(3.31) 
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Our final result for M? +1 _ t \ is 



-is 



3-87rea 2 s f p \t\(M) 



y/2NQ*m* 



(3.32) 



Perfoming the integral over z in ( 3.1 1|) we find 



M, 



(0,0) 



7ia s Qm 2 



47^±ilxG(x, QV4) 



(3.33) 



Now let us discuss M} +1 _ x y Using ( |2.11|J3.2| ) and the second term of ( |3.6| ) we 
can write 



M ( + i,-i) 



is 



7*P 



d 2 rdzr 5 3ef p 
2tt 8^12iVV2 



{M)Q(zzf /2 K 1 (rQVzl)I 1 . (3.34) 



In this case the integration over the transverse momenta of the t— channel gluons 
is different as compared to M? +1 _^ case 



h 



a s S a 



d 2 k 



[2fc 2 (kq) 2 - k 4 q 2 - 2(kq) ; 



f k 2 (k-q) 2 

3(kq) 2 g 2 (z 2 + z 2 ) - 2(kq)g 4 (z 3 + z 3 )} J p ab (k, q) . (3.35) 
Extracting the main, logarithmic, part from the above integral we find 



a s S ab q A [3(z 2 + z 2 ) - 1] f d 2 k 



h 



(3.36) 



q 2 2 

Therefore we can relate M} +1 using the relation ( |3.9| ), to the gluon density. 

oo 

Performing the integral over r, using / Ki(r)r e dr = 384, we have 



ML 



(+i,-i) 2S 7 * P 



dz 



A-3TT 2 ea s f P \t\{M) [3(z 2 + z 



V2NQ e 



zz 



xG(x,Q 2 zz) . (3.37) 



Note that the z integration is different for Mh x _ x \ and M° +1 _ x y The additional 
factor (zz)Q 2 appears in the denominator of (|3.37|) by virtue of two additional 
powers of r in the numerator of eq. ( |3.34| ). As a result the z integration for 
Mh x _ x \ becomes similar to ones for Mh 1>+1 \, ML +1 ^ and Mh 1>0 y Calculating 
the integrals over z ( using (|3.14j ) in eqs. ( |3.11| ) and ( |3.37| )) we find 



Ml +1> _ x) _ \t\(M)2(j + 2) 



Mt 



(0,0) 



Q 2 



7 



M, 



Mt 



(0,0) 



rf + r] 1 



(3.38) 
(3.39) 
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It should be noted that although Mh x ^ is suppressed as compared to 

2 

_i) by the factor ~ ^f, the perturbative part M^ +1 _^ can be dominant 
at high energies since it has a steeper energy dependence (~ xG(x, Q 2 )) as com- 
pared to the soft part M? +1 M} +1 _^ contains also an enhancment factor I/7 
which originates from the more singular z integration. We will show below that 
\Mh x _ X J > M? +1 _ x % at typical for HERA kinematical conditions. 



3.4 Additional remarks on the amplitudes 

We have calculated above the helicity amplitudes for the diffractive vector meson 
electroproduction at s 7 » >> Q 2 >> Aq CD (large Q 2 , small x). Let us discuss now 
the assumptions used and the physical issues related with the helicity amplitudes. 

We assume that perturbative QCD can be applied to describe the process 
at large Q 2 . Although the QCD factorization theorem has been proven 
only for the leading scalar amplitude M(o,o), we extend here the perturbative 
aproach, following (where transverse helicity non-flip amplitude M( + i )+ i) was 
considered), to describe helicity-flip amplitudes. 

The main features of the perturbative QCD aproach are the following. The 
virtual photon splits into the massless qq pair having the total helicity 0. As a 
result, the helicity of the photon coincides with the z projection of the quark 
angular momentum. The helicity states of the quarks do not change during the 
interaction with the proton. Therefore the helicity of the meson is equal to the 
projection of the angular momentum of the outgoing qq pair onto the direction of 
the meson momentum. The helicity flip comes from the change of the projection 
of the qq angular momentum during the interaction. This change originates in 
our aproach from the non-forward kinematics (t 7^ 0). Therefore there is no 
suppression of the helicity-flip amplitudes with energy, they are driven by the 
leading gluon (pomeron) exchange. M( + i i0 ), M(o,+i) and M} +1 _^ have the energy 
dependences which are similar to the energy dependences of the helicity non-flip 
amplitudes. They are proportional (at small t) to xG(x,Q 2 ). 

According to our calculations the helicity single-flip amplitudes are propor- 
tional to the double spin-flip one is ~ \t\. These factors come from the 
expansion of the 7* — > p transition impact factors describing these amplitudes. 
Since the typical transverse separation between the quarks is small, ~ Q^/f > this 
expansion is expected to be valid up to the rather large values of \t\, \t\ < Q 2r y. 



On the other hand, equation ( |3.9|) relating the proton impact factor with the 
gluon density and its implications eqs. ( p.ll| , |3.13j|3.17| , |3yT9| , |3.37|) are valid only at 



very small (vanishing) momentum transfer. The coupling of the two gluon system 
with the proton decreases with the growth of the momentum transfer. We are not 
able to describe the t behaviour of this coupling from the first principles. But we 
see that at small t the structure of the integration over the transverse momenta 
of the t— channel gluons is similar for both the helicity non-flip and the helicity- 
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flip amplitudes (with the exception of the nonperturbative part of the double-flip 
amplitude M? +1 _^). Therefore it is natural to expect an universal t behaviour 
for all helicity amplitudes coming from the coupling of the t— channel gluons 
with the proton. And we believe that this universal t dependence is canceled in 
the ratios ( |3 . 1 5| , |3T2"0"| , |3 . 2 1| , |3~38| ) and, therefore, these results are valid in the broad 



t region extending up to \t\ < Q 2r )- This assumption is based on the observation 
that in all cases we deal with the scattering of a qq pair the transverse size of 
which is much smaller than the size of the proton. 

We will parametrize this universal t behaviour as ~ e -6 '*'/ 2 for the amplitudes 
(~ e~ fe '*' for the cross sections) with the slope b which is of the order of the 
square of the proton size. Due to the large value of this slope, b = 5 ... 6 GeV -2 



according to the HERA data |T0|, the helicity flip amplitudes are peaked as well 
as the helicity non-flip ones at small t (\t\ < to ~ i n spite of the fact that 
they contain the factors J\t\ for single-flip and \t\ for double-flip. 



According to eqs. (|3 . 1 5| , |3 . 2 0| , |3 . 2 1|J3 . 3 9f) for the typical \t\ values , \t\ « 1/6, 



assuming that (M) ~ m p , we have 

1 > a > (3 > \5\ > \rj\ . (3.40) 

Therefore the largest among the amplitudes violating SCHC, the helicity single- 
flip amplitude M( + i j0 ), is smaller than the transverse non-flip amplitude M(+i )+ i). 

Let us estimate the ratios in eq. ( j3.40| ) for the kinematical conditions relevant 
for the HERA experiments. We choose Q 2 = 10 GeV 2 , x = 10~ 3 which corre- 
sponds to W — y/s^p = 100 GeV. We will give the estimates for the simplest 
situation when the momentum transfer is not restricted during the helicity anal- 
ysis, i.e. the experimental sample is not divided into the t bins. In this case we 
can substitute the factors J\t\ and \t\ by their mean values 



(1*1) 



dt\t\e~ m 



1 



(Vl*l> 



du[\t\e~ m 



J dte ~m ~ 2Vb 



We shall use 6 = 6 GeV -2 , (M) = m p in our estimates. For the effective gluon 
anomalous dimension at these values of Q 2 and x we use two values 7 = 0.7 and 
7 = 0.5. These values of 7 are close to the ones presented in Fig. 3 of 0. 
The results of our estimates for 7 = 0.7 (7 = 0.5) are the following: 

a = 0.59(0.73) , (3.41) 

P = 0.12(0.16) , (3.42) 
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5 = 0.056(0.079) , (3.43) 

r] 1 = 0.031(0.041) , (3.44) 

rf = -0.016(-0.015) . (3.45) 

Estimating 77° in eq. ( |3.33| ) we use as = 0.3, as = 0.87. This value of as = 0.87 
used in eq. (|3.29|) reproduces the correct value of the total pp cross section at 
^/Sp~^ = 100 GeV. Combining eqs. ([3.44 , 3~I4"5|) we obtain a very small number for 



the ratio of Mf+i-y to M( 0i0 ) 

77 = 0.015(0.026) . (3.46) 

M( +li o) is approximately 7. . . 8 times smaller than Mmo)- The other helicity flip 
amplitudes are considerably more suppressed. 

We have calculated above the dominant at high energy imaginary parts of 
the helicity amplitudes. We will give here only an argument why the real parts 
can be not too important in the polarization phenomena. The real parts of the 
amplitudes are related to the imaginary ones through the dispersion relations. 
Since the imaginary parts of all helicity amplitudes have according to our con- 
sideration similar energy behaviour we expect that the helicity amplitudes have 
phases which are close to each other. Therefore observable effects related to the 
differences of these phases will be additionally suppressed. 



4 Vector meson decay angular distribution at 
HERA kinematics 

The polarization of the p meson is experimentally accessible through the mea- 
surement of the angular distributions of the decay products. For the relations 
between the helicity amplitudes and the angular distributions we shall use the 



results and standard conventions of ||, see also |TIJ. The definition of the three 
independent angles involves three planes: 1) the electron scattering plane, 2) the 
vector meson production plane (which contains the photon and meson momen- 
tum vectors), 3) the meson decay plane. The orientation of the meson decay 
plane is described by the polar and the azimuthal angles (9) and (0). The third 
angle ($) is the angle between the electron scattereing and the meson production 
planes. 

The polarization parameter of the virtual photon density matrix, 

- l ~y 

e 2 ' 

l-y + \ 

is close to 1 at HERA kinematics. For W = 100 GeV (y = = 1/9) its value 
is e = 0.993. 
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The decay distribution W (cos 6, (ft, contains the parameters, the matrix 
elements rf k , which are known bilinear combinations of the helicity amplitudes 
0. These matrix elements can be determined experimentally by the analysing 
moments of the observed decay angular distribution. 

According to our estimate the helicity flip amplitudes are substantially smaller 
than the helicity non-flip ones. Nevertheless, as we shall show, the largest among 
them M( +lj0 ) leads to a sizable effect. 

The ratio of the longitudinal to the transverse cross section is expressed 
through the ratios of the helicity amplitudes as follows 

R = a L /a T = N L /N T , (4.1) 

where 

N T = a 2 + 1 + rf ■ N L = 1 + 25 2 . (4.2) 

The matrix elements entering the decay angular distribution have the following 
expressions in terms of the ratios of the helicity amplitudes ( assuming that the 
helicity amplitudes are purely imaginary): 

r°* = B(e + f3 2 ) , Re(r%) = B(2e5 + /3a - {3 V )/2 , = B(a V - e5 2 ) ; 
r\ x = Bar, , Re(r\ ) = Bf3( V - a)/2 , = -Bf3 2 , r\_ x = B{a 2 + V 2 )/2 ; 
Im(r 2 w ) = B(3(a + rj)/2 , Jm^) = B(rf - a 2 )/2 ; 

= fy(a - V) , Re(r 5 w ) = -^(2{35 + a- V )/2 
r 5 = -|2/3 , rU = -^5( V - a) ; 
/m(r 6 o) = __B (a + v y 2 , Im{rl_ x ) = -^6(a + V )/2 . (4.3) 

We introduce for short the notation 

B = \/{N T + eN L ) . 

Substituting our estimates for the ratios of the helicity amplitudes derived in 
the previous section we have 

r{$ = 0.74(0.65) , Re(r%) = -0.015(-0.014) , 
r f_ x = 0.0042(0.0082) ; 
r\ x = 0.0065(0.012) , Re(r\ ) = -0.025(-0.036) , 
r i Q = -0.011(-0.016) ,r\_ x = 0.13(0.17) ; 
Im(r 2 ) = 0.027(0.039) , Im(r 2 ^) = -0.13(-0.17) ; 
r 5 n = -0.017(-0.025) , Re(r lQ ) = 0.15(0.15) , 
r^o = 0.12(0.14) , rj_! = 0.017(0.025) ; 
Im(r%) = -0.16(-0.17) , /m(rf_ x ) = -0.0088(-0.013) . (4.4) 
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The matrix elements which have nonzero values in the case of SCHC ( if 
I3,S,ti = 0) are r^, rj_ 1; Im(rf_i), i?e(rf ), Im(rf ). The relations which would 
hold between them in the the case of SHCH are only slightly violated since 6, rj 
are very small: 

-(1 - r{£) - t\_ x = Be6 2 « 0.002(0.004) , (4.5) 

r\_ x + Im{r\_^) = Brj 1 « 2(4) • 10~ 4 , (4.6) 

Re(r 5 10 ) + Im(r 6 10 ) = -^{(35 - rj) » -0.011(-0.017) . (4.7) 

v2 

Therefore these matrix elements should be measured with high precision to see 
the violation of the SCHC. Note also that the value of R calculated using the 



„04 

SCHC relation R = ~ - ' 00 04 exceeds the one calculated using eq. ( |4.1| ) by 

1 — r oo 

3... 4%. 

It is natural that the effects of the violation of SCHC manifest themselves 
more transparent in other matrix elements which would be zero in the case of 
SCHC. The largest among them is Tq « 0.12(0.14). This matrix element has a 
very clear meaning. It is related to the interference between the two amplitudes 
describing the two possibilities to produce the longitudinally polarized vector 
meson. These are the dominant at large Q 2 amplitude M(o,o) and the helicity 
single-flip amplitude M(+i f o). 

Since this interference exists on the level of the production of the vector 
meson and does not depend on the kinematical variables describing the meson 
decay (9,<f>), the correspondig effect survives without the loose of the analysing 
power after the integration of the angular distribution over the angles (6, (f>). The 
resulting distribution over the relative angle between the electron scattering and 
the meson production planes is 



W{$) 



l + v /2e(l + e) cos $(^ + 2^) 



(4.f 



We skipped in the above equation the term ~ cos 2$ which is proportional to the 
small matrix elements Tq , r\ v Substituting our estimates for Tqq, we obtain a 
substantial deviation of the $ distribution from the flat one 



W{$) = [1 + 0.18(0.19) cos$] . (4.9) 

The other matrix element that could be potentially large is i?e(r°o)- Since 
it contains the term ~ 2eS which is linear in 6, this matrix element has a large 
sensitivity to the second single-flip amplitude M( 0) +i). But its value turns out 
to be small due to a large cancelation between the terms 2e<5 and j3a. This 
cancelation is related with the opposite signs of M( +lj0 ) and M( 0j+ i). Note that in 
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the case of positive sign of 5 this matrix element would be estimated as Re(r 
0.066(0.086). 



on 
10 ) 



The matrix elements (O) parametrize the angular distribution for the un- 
polarized initial positron. A few additional matrix elements can be measured if 
the initial positron is longitudinally polarized. But these matrix elements do not 
exhibit an advantage in sensitivity to the violation of SCHC as compared to the 
unpolarized ones. 



5 Conclusions 

We have considered in perturbative QCD the polarization effects in diffractive p 
meson electroproduction. We assume that perturbative QCD with the account 
of the important effect of the gluon scale behaviour is applicable to all helicity 
amplitudes. We estimate the helicity amplitudes in the approximation of the 
constant effective anomalous dimension of the gluon. Our results are summarized 
intheeqs. (|37T5|JCT| , F21| , |338|j339| ). 



The equations derived here can be applied also to other light vector mesons 
lu, (j) with corresponding changes in the coupling constants and parameters (M) 
describing the meson wave functions. 

The perturbative QCD leads to a very definite qualitative picture for the 
violation of SCHC at high Q 2 . The largest among the helicity-flip amplitudes is 
M( +10 ). At typical values of t this amplitude is smaller than the transverse helicity 

non-flip amplitude M( +1+1 ), M( +1>0 ) ~ -^^M( + i )+ i) . The other independent 

single-flip amplitude M( ,+i) is suppressed compared to M( + i )0 ) by the factor -jj 1 - 

The double-flip amplitude M( +1 _!) consists of the two parts. The perturbative 
contribution M} +1 ^ is larger at small x than the nonperturbative one. But it 

is suppressed compared to M( ,+i) by the additional factor ~ Therefore 

at very high Q 2 we have M( +10 ) >> |M( 0) +i)| >> M( +1 _!). For the kinematical 
region typical for the HERA experiments we find that |M( 0j +i) | is about two times 
smaller than M( +1)0 ), and M( +1 _x) is about 10 times smaller than M( +10 ). 

This hierarchy between the helicity amplitudes leads to the peculiar predic- 
tions for the parameters of the meson decay angular distribution. We predict 
that the only one parameter among that vanishing in the case of SCHC, Tq , 
deviates substantially from zero. The relations between the parameters, which 
are nonzero in the case of SCHC, are only slightly violated. It would be very 
interesting to confront these predictions with the data. 

Let us note that the parameter Tq ~ M^im is sensitive to the meson wave 
function. Since single spin-flip transitions are proportional to the factor vanishing 
if z = z, see ( f2.22|JT8|) , the helicity amplitude M( +10 ) whould be zero if the wave 



22 



function of the meson is a non-relativistic one ~ S(z — 1/2). It is a consequence 
of perturbative QCD that the fluctuation of the light meson into the pair of 
current quarks is described by the broad wave function. On the other hand in the 
nonperturbative models a meson is often considered as a weakly bounded system 
of constituent quarks having a mass m q ~ ttim/^ and, therefore, described by a 
function close to 8{z — 1/2). Therefore the large value of 7q is a characteristic 
prediction of perturbative QCD. 

Note that our work is only a first step in the investigation of the helicity-flip 
amplitudes. More detailed numerical calculations can be done using the equations 
derived in this paper. It is possible to do this without the approximation of 
the constant gluon anomalous dimension and to investigate in more details the 
dependence of the amplitudes on the meson wave functions. Also the real parts 
of the helicity amplitudes have to be considered. 
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